6 years of harvest, 2020 had variable rate seeding
plot.gps.grids(corn.2013.dat,plots=0,residuals=TRUE)
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
plot.gps.grids(soybean.2014.dat,plots=0,residuals=TRUE)
plot.gps.grids(corn.2015.dat,plots=0,residuals=TRUE)
plot.gps.grids(soybean.2016.dat,plots=0,residuals=TRUE)
plot.gps.grids(corn.2017.dat,plots=0,residuals=TRUE)
plot.gps.grids(soybean.2018.dat,plots=0,residuals=TRUE)
corn.2013.var <- variogram(Yield~1,
locations=~Easting+Northing,
alpha=c(0,45,90),
data=corn.2013.dat)
plot(corn.2013.var)
soybean.2014.var <- variogram(Yield~1,
locations=~Easting+Northing,
alpha=c(0,45,90),
data=soybean.2014.dat)
plot(soybean.2014.var)
corn.2015.var <- variogram(Yield~1,
locations=~Easting+Northing,
alpha=c(0,45,90),
data=corn.2015.dat)
plot(corn.2015.var)
soybean.2016.var <- variogram(Yield~1,
locations=~Easting+Northing,
alpha=c(0,45,90),
data=soybean.2016.dat)
plot(soybean.2016.var)
corn.2017.var <- variogram(Yield~1,
locations=~Easting+Northing,
alpha=c(0,45,90),
data=corn.2017.dat)
plot(corn.2017.var)
soybean.2018.var <- variogram(Yield~1,
locations=~Easting+Northing,
alpha=c(0,45,90),
data=soybean.2018.dat)
plot(soybean.2018.var)